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Abstract 

The Dirac equation for an electron in an external electromagnetic field can be regarded 
as a singular set of linear equations for the vector potential. Radford's method of al- 
gebraically solving for the vector potential is reviewed, with attention to the additional 
constraints arising from non-maximality of the rank. The extension of the method to gen- 
eral spacetimes is illustrated by examples in diverse dimensions with both c- and a-number 
wavef unctions. 
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1 Introduction 



The Maxwell-Dirac equations are the coupled nonhnear partial differential equations 
which describe a classical electron interacting with an electromagnetic field. They are 
also the equations from which quantum electrodynamics is derived. Since the mathemat- 
ical foundations of the latter remain unclear, the Maxwell-Dirac equations continue to be 
of interest [|I|, ^ ^. Recently Radford^] handled the Maxwell-Dirac equations firstly by 
solving the Dirac equation for the electromagnetic potential in terms of the wavefunction 
and its derivatives, and then substituting this solution in the Maxwell equations. This 
approach subsequently led to some physically interesting results H, (for a review see 

§)• 

Despite the viability and potential importance of Radford's algebraic solution, at least 
for the treatment of the equations of classical electrodynamics for electrons and photons, 
it appears that the method has not before appeared in this context. One analysis which 
reached negative conclusions about the approach, and which may have engendered the lack 
of attention to it in the literature, is that of EliezerQ. In that paper, it was noted that the 
determinant of the matrix of coefficients of the vector potential in the Dirac equation 
actually vanishes, and that therefore a umique algebraic inversion was not possible. The 
aim of the present paper is to reconcile 0] and @ , and to emphasise the legitimacy of the 
algebraic Ansatz, despite the negative conclusions of 0, by a careful analysis of the nature 
of the Dirac equation regarded as a linear system for A^. The main result is that the 
Dirac equation is indeed invertible if a real solution for the vector potential is required, 
and moreover that the treatment entrains an additional set of polynomial constraints on 
the wavefunction and its partial derivatives which must be carried forward in any further 
analysis. In §2 below, the abstract formalism is developed, and (for the 4 dimensional 
case) it is shown how the explicit manipulations, which rely on the structure of the Dirac 
algebra to derive the solution for the vector potential and the additional constraints, 
conform to the general setting (it is also pointed out that the solution can be regarded 
as including the mass, or more generally a Lorentz scalar potential, as a fifth unknown). 
This is done both in Lorentz-covariant Dirac spinor notation, and in van der Waerden 
2-spinor notation. In §3, the case of arbitrary (flat) spacetimes with signature [t, s) is 
taken up. Known results on the structure of the Dirac algebra (formally, the Clifford 
algebra C{t,s)) in these cases are used to give an enumeration of constraints which are 
quadratic in the wavefunction and derivatives (in addition to current and partial axial 
current conservation, which hold in all cases). The four- dimensional results are recovered, 
and generalised to the case of a-number as well as c-number wavef unctions. A major 
outcome is a tabulation (table ^ of such constraints as to fermion wavefunction statistics 
and metric signature in diverse dimensions. Concluding remarks and prospects for further 
development of the work are given in §4 below. 



2 The 4 dimensional Dirac equation 



8x4 real system 

The Dirac equation for a fermion of charge q described by the spinor wavefunction il) in 
the presence of an external electromagnetic potential may be written^ 

q-i^A^ = - mi)). (1) 

Following Eliezer[^], we write this as a matrix equation for A^] 

M^A, = (2) 

where = Yj^fs for M e M^iC) , M : — ^ and A, Z e 

In it was noted that M had rank 3 and determinant zero, with a rank 1 right null 
space, and therefore could not be inverted to obtain a unique solution for the potential. 
Yet RadfordQ did just this, albeit in the bispinor representation. That work exploited 
the fact that is realQ, which was not used in 0. The point is that despite the zero 
determinant, (|^) can be inverted if we know that A^ is real, and providing that the 
intersection of the right null space of M with (as a subspace of ) is trivial. Even 
though the columns of M are not linearly independent as a vector space over C, they are 
in general linearly independent as a vector space over M. 

We may break (0) into real and imaginary parts, yielding a system of 8 real equations 
in 4 real unknowns, schematically AiA = Z, where 

i(M + M*) \ __ / + Z 



^-[l(M-M*))^ ^-[liZ-Z*))- 

Ai is not square; no determinant is defined; yet there are other tests for linear inde- 
pendence and invertibility||10||. To invert a system of m equations in n unknowns, with 



m > n, of the form (^, then we seek an n x m matrix Q such that QA4 = 1, the unit 
n X n matrix, and then A = QZ. If such a solution exists, the rank of is n (= 4). 

Note that the multiplication of a row of Q with a column of M. is actually a Euclidean 
real inner product. If the columns of M. are understood to be spinors over an 8 dimensional 
real basis, we can accept the same interpretation for the rows of Q. The existence of an 8 
dimensional real basis thus supplies us with a definition of a real inner product between 
spinors, for which we will use the notation {(p • ip). It is easy to verify that (</> • ip) is 
actually equal to the real part of the standard complex inner product: 

{(f).x) =^e<(f),x>=^{(f>^X + i^^x)- (4) 



^ In this section standard Cartesian coordinates x^^ // — 0,1,2,3 for four dimensional Minkowski 
space with (1,3) metric ('/7/iL')/i,i^— 0.1,2,3 = diag(+, — , — , — ) are introduced. AfRces for Dirac spinors are 
introduced as ct = 1,2, 3, 4, while the Dirac matrices (generators of the Chfford algebra C(l, 3) in the 
standard basis) satisfy {7/^,7;^} — 2r]^^; for conventions sec In §3 below, the notation is generalised 
to dimension d = t + s with metric signature {t, s). 

^ The most telling use of the reality of Afj_ is implicit in the bispinor representation, where half the 
equations are conjugated, taking to be real. This results in systems of equations where the matrix 
operating on can have non-zero determinant |7|. This transformation cannot be performed with a 
matrix transformation on C*, and does not preserve the determinant. 



The system of m equations in n unknowns entails that the right hand side of the 
equation should satisfy m — n (8 — 4 = 4) additional consistency conditions, arising from 
the fact that Z must fall in the column space of Ai. To find these consistency conditions, 
we seek a further m — n linearly independent spinor rows x that have zero real inner 
product with the columns of A4 {x span the left null space of A^). The consistency 
conditions may then be written {x • Z) = 0. Q is not unique, for any linear combination 
of the rows x can be added without changing its effect on M.. 

It is not necessary to work explicity in 8 real components: regardless of which basis 
we use, the columns of M. and the rows of Q are just spinors in a vector space over M, and 
the matrix multiplication is just the calculation of inner products using (|^). All that we 
require for the inversion is to find spinors (py where • 'y'^ip) oc 6^. For the consistency 
conditions, we require 4 linearly independent spinors x such that {x • 'y'^ip) = 0. As will 
now be shown, the structure of the Dirac algebra indeed admits such rows, (pu and x- 



Inversion 

Let = '~i^'~iu'ip. Then 



We use the hermiticity of 7° and 7^7^: 

and likewise for 7,^. Then 

(0, . ri^) = i(V;t70(7,7'^ + 7/^7,)^;) = 

where the Dirac conjugate ip is defined in the usual wayQ as ^/^ = ip'^'~^^ . 

Applying the real inner product with the same rows to the right-hand side of the Dirac 
equation (|I]), gives explicitly 

The last term is identified with the current j'^ = ip'^'^ip-, so we can write the solution for 
the vector potential 

^ ^ 1 ^(V^7m# - ^ 7/.V^) - 2m 



General properties and nomenclature for the Dirac algebra in generic dimensions and spacetime 
signatures (including four dimensional Minkowski space) are given in §3 below. 



Consistency Conditions 



It is possible to show by the use of (H) two separate sufficient conditions for spinors x to 
have zero real inner product with 7^7/;: 

1. X = r-?/;, where F is a matrix in the Dirac algebra such that r''"7'^ is an^zhermitean; 
Alternatively, 

2. X = Tip* where T is a matrix in the Dirac algebra such that T^7^ is antisymmetric. 

As an example of (1), take x = i'y^'ip- Both the left hand side as well as the mass term of 
(|l]) vanish (as Z7° is itself antihermitean) , leaving 

This is the normal current conservation equation, 

d-j=d^f = d^ilj^^''Y^ = 0. (6) 

Also satisfying condition (1), take x = using the antihermiticity of i7°757'^. The 

hermiticity of ^757° now ensures that the mass term survives, and manipulations on the 
right hand side of (|lD lead in a similar way to the equation for partial conservation of 
axial current jg = ip'j^'j'^ip as the second consistency condition: 

(9 ■ J 5 + 2imip'~^^%lj = 0. (7) 

As an example of the sufficient condition (2), take x = InCip* and x = il^Cip* 
respectively, where C is the charge conjugation matrix. We evaluate these inner products 
with '-y^'ip using the hermitean conjugate {'^^Cip*)'^ = ip^i'JdCy = —ip^Cj^. These inner 
products are the real and imaginary parts of —ip^Cj^'j^ip respectively, which is zero by 
the antisymmetry of 6*7^7'^. Applying the same row operation(s) to the right hand side 
gives = ip^C'y^li'y^dfj^ip — mip) or 

iP'C'y^^iP = (8) 

(after using the antisymmetry of C75 itself to eliminate the mass term) yielding one 
complex condition, or two real conditions on the spinor. This is the result previously 
reported by EliezerQ (who attributed to Dirac the antisymmetry argument using C75 = 
a^Qiz in the standard representation). The consistency conditions (^, (|^), and (^ are 
equivalent to Radford's[§| 'reality' conditions. 

Alternative inversion 

As mentioned above, the choice of nonsingular matrix inverting Ai, and consequently the 
form of the final expression for A, is not unique. As an alternative choose (pi, = i'y^'y^'ju'ip- 
We then find by a similar working to (0), using the anticommuting property of 75 with 
that _ 



Applying the inner product with the same rows to the right hand side of the Dirac equation 
(HD, in this case the mass term vanishes, yielding an alternative solution for the vector 
potential: 

That (|) and are indeed equivalent, and equivalent to follows from the use of Fierz 
identities together with the auxiliary constraints (see §3 below). 

8x5 real system 

The inversion (y) does not contain any mass term. However, note that the pseudoscalar 
consistency condition can be written 

m = = . (10) 

The similarity between (^ and ([10|) suggests that the original system could have been 
considered as 8 real equations in 5 unknowns, qAq, . . . , qA^, and m (or more generally a 
Lorentz scalar potential). In this system, (^ and (0) provide an inversion, while (^ and 
the real and imaginary parts of (H) provide the 3 consistency conditions. 

2-spinor analysis 

Radford and Booth and Radford |§ used van der Waerden notation in order to derive a 
complex form of the vector potential subject to additional reality conditions. Here the 2 
spinor version is reached via the Weyl representation of the Dirac algebra (see for example 
[pnp, wherein 



V J \U J \ Un 

A generic matrix F in the Dirac algebra has matrix elements 



in particular. 



where0 



a 



^f^ab Q 

{crno<,<3 = {cr^ cr), (a^)o<M<3 = {cr' , -(t). (11) 



''The Pauli matrices are 



£ab 
E^'' 



The definitions (0) are consistent with 

Cap 

together with e = ia'^, that is, component-wise 
Starting then from 



C 



'ab 



1 
-1 



ab 



—e 



ab' 



and transcribing to 2-spinor form, the Dirac equation reads directly 



where0 



A 



ab 



a 



M A 



A 



ab 



a^^^A^. 



Finally taking complex conjugates. 



u ■ 


qA^'^-v, = 


-id'^'^v^ + niu" 


V ■ 


qA^\b = 


id^^Ub + mv\ 


u ■ 


qA'% = 


id'"^Vd + mu", 


V ■ 


qA''^^ = 





Thus by taking combinations of the form a{V"'v — U"'u ), f3{U u"" — V v"") and using 
UdV^ — VdU^ = 5d'{ucV^) the vector potential can be isolated, with general 'solution' 



A 



cd 



12 



with arbitrary parameters a, (3. As emphasised by Radford0], all these forms are equiva- 
lent, subject to the hermiticity conditions satisfied by the potential itself. The latter can 
be imposed via the two-spinor projections of A"*, namely 



{A^'UaUb) 
{A'^'VaVb) 
{A^^^UaVb) 



{A^'UaUb) 
{A'^'VaVb) 
{A^'VaUb) 



Substitution of these conditions into a suitable form of 
0, leads to 



(13) 

for example with a = 1, (3 = 



t^ab 



0, 



uXbd.v' 



2im{v°'Ua - UaV") 

Vaa^^'d^Ub. 



(14) 



^Hermiticity, and raising and lowering of indices are entailed in the relations cr^^ = £afc£(jj,f^ i 



Referring ( ]T^ ) to a fixed tetrad basis, UaUb, VaVb, UaVh and VaUh via appropriate contrac- 
tions and using (0) tlien shows directly that all forms are equivalent to the manifestly 
hermitean version with a = (3 = 1. As expected from the general analysis in the previous 
subsection, ([T^ ) with a = [3 = 1 agrees with the previous result (||) expressed in the Weyl 
basis, and ([T^) are equivalent to current and partial axial current conservation^ and the 
additional complex pseudoscalar identity (|]) satisfied by the Dirac wavefunction. 



3 Higher dimensional extensions 

In the previous section it was shown that the 4 dimensional Dirac equation can be regarded 
as a singular set of real linear equations for the vector potential A^. Gaussian elimination 
in this 8x4 real system then leads to a solution for and also implies a set of four 
additional linearly independent constraints, linear in 9^^, which can be identified in this 
case with bilinear identities proposed by Radford and Booth and Radford using van 
der Waerden 2-spinor notation. Either approach ultimately derives from the structure of 
the Dirac algebra and the symmetry properties of the 7 matrices. The same analysis is ex- 
tended here in Dirac notation to higher dimensional cases and different metric signatures 
(in fiat spacetime). Also, an important distinction to make is that between c-number 
and a-number (or Grassmann-valued) Dirac spinors. For the latter, a Gaussian elimi- 
nation argument requires a formal treatment of linear algebra over Grassmann-extended 



ground fields |T2|; for present purposes it suffices to assume that the count of solutions 
and constraints goes through. 

The explicit construction of the Dirac algebra in higher dimensions has been reviewed 



by Tanii|jT3|; see also |ll4| , |l5t The Dirac matrices satisfy 

where the spacetime metric signature is taken as (t, s) for even dimensions t+s = d, indices 
being labelled /i, = 0, 0', 0", . . . ; 1, 2, . . . , s, with chirality determined by the equivalent 
of 75, 7 = 7o7o, . . . 7i72 • • • 7d- 

There are three involutive automorphisms associated respectively with complex con- 
jugation, transposition and hermitean conjugation under which the 7 matrices acting on 
complex spinors ipa undergo similarity transformations^, 

(16) 
(17) 
(18) 

where Sa, Sb, Sc are sign factors depending on the spacetime: 

Sa = 5b5c ; 



7 ^ 


= Sa 


A , 


Ifi pi 


{A- 




7 ^ 


= 5b 


-Da 


■y* 1^' 




J/3' 


7 " 


= 5c 


Caa' 


Ifi pi 


{c- 





^ Take the sum and difference of (|l^) and ( 14 o) 

^The notation "yii \ 'yjf^ is used for the expressions on the right hand sides of (fsl) (without the 



sign factors), extended to arbitrary elements of the Dirac algebra T (see below). Note that in the previous 
section the index conventions for the A and C matrices differ from that used here in the general case. 



Table 1: The values of as a function of 5b = ±1 and s — t (mod 8) (from [O; x 
indicates that no representation with the specified signs exists). 



s-t 





1 


2 


3 


4 


5 


6 


7 


e+ 


+ 1 


X 


-1 


-1 


-1 


X 


+1 


+1 


e_ 


+ 1 


+1 


+1 


X 


-1 


-1 


-1 


X 



5b = ±1; 



6r: = 6u(-lY+\ 



A, B and C are related through the definitions of Dirac and charge conjugation for spinors, 



SO that 

The fundamental identity B*B = £^11 determines the existence of Majorana spinors 
for metrics in which = +1, and 5b = —1- Dimensions for which this is possible can 
be read from table which gives the values of es = \/2cos[^7r(s — t — 5b)]- Finally, 
the implementation of parity in the Dirac algebra is determined by the symmetry of C, 
C^ = 6cC, ec = {5B)\-lf^'^'''^eB. 

Consider now the Dirac equation for ip and conjugate forms: 

qJ^ip = — rnjip 
q^-"^" = (—■ — 5Bm)ip'' 

= '^{i^ - 5cm) (19) 

where the sign factors follow from the explicit definitions in (|T^). Write (|T^ as ^^'^ 
and '^^ respectively. Then, as in the four dimensional case, the solution for the potential 
follows by taking combinations ipj^"^ + ^T^"?/" which force the isolation of through an 
anticommutator of 7 matrices: 



2q ipip 



where sa = |(1 + 5a) = ^(1 - (-1)*) (compare (||)). 

The identification of additional identities satisfied by the Dirac wavefunction given 



(19), (|20|) amounts to determining the structure of a certain ideal in the free algebra of 
rational expressions in ip, ip* and partial derivatives d^ip, d^ip*, for appropriately smooth 
wavefunctions. Such expressions do not necessarily form linear representations of the 



spacetime Lorentz symmetry group SO(t,s), and the distinguished role played by the 
constraints following from Gaussian elimination is not clear. Here, the quadratic case is 
analysed, by analogy with the four-dimensional case, and in relation to the counting sug- 
gested by dimensional considerations. Polynomials in ip, ip and derivatives do decompose 
with respect to the Lorentz algebra. Moreover, there is a natural bi-grading by degree: 
in the quadratic case, this is simply by fermion number F. Thus the cases ip ■ ip, with 
F = 0, and tp ■ tp, or equivalently ip'^ ■ ip, with F = 2, can be considered (as the latter will 
necessarily be complex, constraints in ip ■ ip do not require separate treatment). A case 
by case analysis follows. 

F = ±2 : 

As is well known, products of two spinors admit a decomposition into antisymmetric 
tensor representations of the Lorentz group [0. Noting that 

= eBiC-'vr^^^^f,, (21) 

it is clear that an explicit decomposition is provided by the linear basis for the 
Dirac algebra comprising the antisymmetric p-fold products of 7 matrices [[13|, |18[ 
{7m, Ifiu, 7mim2...mp' ■ ■ ■}, P = ^, ■ ■ ■ ,d. Furthermore, for appropriate p, depending 
on the spacetime signature, dimension and statistics of the spinor wavef unctions, 
(plD is identically zero, leading after contraction with to a sequence of differen- 
tial identities in the Dirac wavefunction, after implementation of (|1^, ^). If the 



statistics of the Dirac wavefunction is specified as ipa'4'i3 = fi'pi'a-, then using 
and 

^(c) - AP+i^ _ xP+i|'_i^|5(p+i)p^ (22) 

we require 

ec/5?.+'(-l)^(^+^)^ = -I. (23) 



Note also I IS 



~ ~(~-'-)^7/ii/J2---/ip7/i + (p4.i)7M/iiAt2---Atp- (24) 



Theorem 1 For spacetime dimension d = Ak, a single Lorentz scalar or pseu- 
doscalar bilinear complex differential constraint of the form 

^4ip = or ^^4ip = (25) 

exists, whenever ec^cf = Tl; respectively. For spacetimes of dimension d = Ak + 2, 
both scalar and pseudoscalar complex constraints hold ifecScf = — 1; otherwise 
neither holds. In general a sequence of Fierz type identities 

4tp = pA[f,^^^jf,^,„f,^]ip, (26) 

holds, where 5p = ^{1 - ec5^c^^ fi-l)^P^P+^^) . Note that are regarded as 

differential conditions on ip via the substitutions (17^), p 



Proof 

Clearly (0), together with the interchange sign factors implies (^) after contraction 
with A^. After substituting (^) and rationalising, the identity is therefore quartic, 
and generically of Fierz type. The only bilinear cases occur when p = (in which 
case the condition ec^cf = ~1 refers to the symmetry of ip'^jjiip), or for p = d (for 
which 7^ and 7 anticommute, and the charge conjugation of 7 determines whether 
the pseudoscalar identity holds). 

□ 

F = : 

Theorem 2 The real type Fierz type identities 

'^P4'l^^w2...^^pi^ = (-l)^V^7/.iM2.../.p4^ + '2pA[^^JP^^^2■■■^^p]^ (27) 

hold, where the appropriate parts of ( [7^ ^ are to be used on the left and right-hand 
sides, respectively. Current conservation and partial conservation of axial current 
hold for all d: 

d^^j^^ = 0, 
d^'^jj^ip + 2imsA(ipJ^P) = 0. (28) 

Proof 

(pTl) follows directly by contraction of (^a) with 
anticommutativity of 7^ and 7 (See (|2^)). 

□ 

Explicit enumeration of these results is complicated because of the multitude of subcases 
involved. In table |^ the counting and nature of the identities is illustrated in diverse 
dimensions, and for given fermion statistics and metric signatures for the Minkowski, 
conformal and Euclidean spacetimes. The explicit expressions for the sign factors in 
the representations of the Dirac algebra available imply the expressions given for 6c€c, 
namely Eb, SbSb and {—1)2'^6b£b for t = 1, t = 2 and t = d respectively. The entries 
within each metric signature class then indicate for which type of fermion wavefunction 
statistics (c-number, / = +1, or a- number, / = — 1 respectively) the indicated scalar 5* or 
pseudoscalar P complex identity exists. Where = +1 exists, bracketed entries indicate 
the choice 6b = —1 or 6b = +1, consistent with the availability either of Majorana (M) or 
pseudoMajorana (pM) spinors respectively, in that dimension and spacetime signature. 

In table 0, the entry for c-number wavefunctions in four dimensional Minkowski space 
corresponding to a single complex bilinear pseudoscalar identity is the original Dirac 
equation case. For a-number fermions, there exists the corresponding scalar equivalent 
(reflecting the properties of 7^ and 7^75 in the Dirac algebra). Taking into account the 
two real constraints (current and partial axial current conservation), there is in either case 
a total of four real conditions, in accord with the count needed to accompany the four 
dimensional vector potential in the 8x4 linear system (see §2). 



28) uses (|T9|) together with 



The counting of constraints in d dimensions must be examined carefully in relation 
to Gaussian elimination in the corresponding 22^'+^ x d real linear system. For example 
in two-dimensional Minkowski space there are no complex (F = ±2) bilinear constraints 
for c-number wavefunctions, and thus only two real (F = 0) bilinear constraints, whereas 
for a-number wavefunctions there are two real plus two complex identities. The apparent 
under- or over- determination of the system (at the quadratic level) must be reconciled 
with the remaining Fierz-type identities (of quartic type). Two dimensions is a special 
case because of the abelian nature of the Lorentz group, but in other cases, the Lorentz 
decomposition of the remaining higher order constraints also bears on the counting. For 
example, in d = 6 (16 x 6 real system) there are again either no or two complex bilinear 
constraints (see table 0). In the latter further 16 — 6 — 2 — 4 = 4 independent 

real conditions are needed. These may represent 4 additional scalar conditions at higher 
order, or a real 6-vector which satisfies two additional conditions equivalent to two of the 
scalar and pseudoscalar conditions. Similar considerations apply to the d = 8 and d = 10 
cases. 

4 Conclusions 

In this work the linearity of the Dirac equation has been exploited as a vehicle to obtain 
an algebraic solution for the vector potential, and previous discussions in the literature 
(0 and [^) have been reconciled. In addition to its role in the problem of obtaining 
(classical) solutions of the full nonlinear Maxwell-Dirac equations, the algebraic method 
has potential application to the nonabelian case, and to Duffin-Kemmer algebras rather 
than Clifford algebras. Modifications to electrodynamics such as the Born-Infeld theory, 
and indeed the nonrelativistic limit of the Dirac equation itself, may also be amenable to 
further study by the method. 

Table 2: Enumeration of the type of complex bilinear identity (scalar, S or pseudoscalar, P) 
admitted by the Dirac wavefunction of the indicated statistics, for various dimensions d and 
metrics of Minkowski (1, d — 1), conformal (2, d — 2) and Euclidean {d, 0) signature (see text for 
details). 
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